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Abstract 

The Hartle-Hawking state of N = 4 SYM at strong coupling and large JVona fixed 
black hole background has two proposed gravitational duals: a black funnel or a black 
droplet. We construct the black funnel solutions that are dual to the Hartle-Hawking 
state on a Schwarzschild black hole and on a class of three-dimensional asymptotically 
flat black hole backgrounds. We compute their holographic stress tensor and argue for 
the stability of these solutions. 



1 Introduction 



A pivotal moment in the study of quantum fields on black hole backgrounds was the realiza- 
tion that black holes will evaporate by emission of Hawking radiation [1] . The fundamental 
issues associated with this discovery have played a central role in the quest for a quantum 
theory of gravity. However, much of this understanding comes from the study of free field 
theories. When the fields are strongly interacting, not much is known beyond a formal 
definition of a Hartle-Hawking state from a Euclidean path integral. 

More recently, the AdS / CFT correspondence [21 [3], H] has provided a means of studying 
certain strongly coupled field theories. States of a large N, strongly coupled field theory in 
a background Ai are mapped to solutions of a theory of gravity in one higher dimension 
with conformal boundary Ai. Typically, the field theory background is Minkowski space. 
In order to understand strongly coupled Hawking radiation, the authors in [HI El E] instead 
applied the correspondence to a black hole background. 

What then, is the gravity dual to a Hartle-Hawking state? From the perspective of the 
field theory, the stress tensor must approach that of a thermal fluid far from the boundary 
black hole. Then in the bulk, the solution should approach that of a planar black hole 
(also known as a black brane). Furthermore, there must be a horizon near the conformal 
boundary that matches up with the boundary black hole horizon. Therefore, the authors 
of [S] conjecture the existence of two families of solutions: the two horizons either join up 
to form a connected horizon, or they are disconnected. In the former, the solution forms a 
black funnel, while in the latter, the solution is a black droplet suspended over a (deformed) 
black brane, see Fig. [TJ Because of the connected horizon, the black funnels are expected 
to be dual to a situation where the boundary black hole couples to the field theory plasma. 
By contrast, the black droplet solutions describe a boundary black hole that does not easily 
exchange heat or other quantities with the plasma. 

A phase transition between these two phases was also conjectured in [5]. If the ratio 
between the size of the boundary black hole and the temperature is small R/T <C 1, then 
the black funnel may have a long narrow neck, and might become unstable to the Gregory- 
Laflamme instability [Sj. It is therefore natural to suggest that black droplets dominate the 
ensemble for R/T <C 1. Conversely, if R/T ^> 1, the droplet horizon might extend deep 
into the bulk until it nearly touches the planar horizon. It may become favourable for these 
horizons to merge, and so the funnel phase might be preferred at R/T ^> 1. 

Various black funnel and black droplet solutions have been constructed from the AdS 
C-metric [HI EJ |9]. These solutions, however, do not have a boundary black hole that is 
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(a) Black Droplet (b) Black Funnel 



Figure 1: Sketches for black droplets and funnels. 

asymptotically flat. A black droplet dual to Schwarzschild was also found in [10J. There, the 
solution does not have a planar black hole, and is therefore dual to an Unruh state^J 

In this paper, we numerically construct the first black funnel solutions with a boundary 
metric that is asymptotically flat. We work in five dimensions where the bulk metric is dual 
to Schwarzschild, and also in four dimensions where we vary the size of the boundary black 
hole. 

We also compute the boundary stress tensor for these solutions. However, in five bulk 
dimensions there will be a conformal anomaly unless the boundary metric is Ricci flat. This 
will generically introduce logarithms in the expression for the stress tensor, which pose a 
difficultly for our numerical method. The four- dimensional case has no conformal anomaly, 
even for non-Ricci flat solutions. It is chiefly for this reason that we do not vary the size of 
the d = 5 black funnels. 

This paper is structured as follows. In the next section, we describe the construction 
of these black funnel solutions and our numerical method. In the following sections, we 
compute the stress energy tensor and study these solutions by embedding them isometrically 
in hyperbolic space. We then finish with a few concluding remarks. 

^^Note that a Boulware state would require a minimum energy configuration, which would be similar to 
the solution found in jlOj . but with an extremal horizon. 
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2 Constructing Black Funnels 



2.1 The DeTurck Method 

Our black funnels require only Einstein gravity and a cosmo logical constant. The action in 
d bulk spacetime dimensions is 

S = Jd d xV=9(R-2A), A = J d - 1 )(d-2) (21) 

where L is the AdS length scale. Rather than solving the Einstein equation directly, we will 
solve the Einstein-DeTurck equation (also known as the harmonic Einstein equation): 

2A 

where £ M = g a/3 (r^ — F£*g) , and f ^ is the Levi-Civita connection of some reference metric 
g. This reference metric is chosen to have the same asymptotic and horizon structure as our 



metric g. For the metrics we are concerned with, (2.2) is elliptic [TT 



In general, solutions to (2.2) will not be solutions of the Einstein equation unless £ = 0. 
In the cases considered here, it has already been shown [TU] that solutions with £ ^ 0, called 
Ricci solitons, do not exist. Therefore, we can use £ to monitor our numerical error; we check 
that £ = to machine precision. This method also has the advantage that a gauge choice 
need not be made a priori. Solving for the equations will also choose the gauge £ = 0. 



2.2 Integration Domain 

The black funnel naturally has a triangular integration domain with three boundaries: a 
horizon, the planar black hole metric, and the conformal boundary. While it is possible to 
solve PDEs in triangular domains, the task is much more cumbersome, and it is more difficult 
to achieve the same accuracy as that of rectangular domains. However, with an appropriate 
choice of conformal frame, this domain can be expanded to a square. 

By careful examination of the region where the bulk horizon meets the boundary horizon, 
one can show that the metric approaches that of a hyperbolic black hok^j In fact, this has 
to be the case for any horizon that reaches the conformal boundary. Close to the conformal 
boundary, and up to terms that only change the holographic stress energy tensor, the ge- 

2 We thank Donald Marolf for this realization. 
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Figure 2: Implicit coordinate transformation from a triangular to a rectangular domain integration 
domain. The point where the horizon meets the boundary is blown into a hyperslice 



where the black funnel line element approaches Eq. (2.3). In this diagram, H represents 
the hyperbolic black hole, d the conformal boundary, P the planar black hole and H the 
bulk horizon. 



ometry is hyperbolic. Furthermore, the most general cohomogeneity one line element that 
is static and manifestly exhibits hyperbolic symmetry for each hyperslice of constant time 
and radius is that of a hyperbolic black hole. A specific hyperbolic black hole metric is then 
fixed by a choice of the bulk horizon temperature and the boundary horizon temperature. 
In this manuscript, we are only considering bulk horizons with the same temperature as the 
corresponding boundary black hole. The only such black hole metric compatible with the 
desired symmetries is given by the zero energy hyperbolic black hole: 



ds. 



L 



2 



Z 2 



dz 2 

-(1 - z 2 )dt 2 + - — — + drj 2 + sinh 2 77 dtt 2 d _ 3 



(2.3) 



Our solutions must approach the r\ — > 00 limit of (2.3). 

Therefore, the aim is to find a solution that has a horizon, approaches the planar black 
hole and the hyperbolic black hole on either ends of this horizon, and has a boundary that 
is conformal to the background metric of the field theory. The above set of considerations 
changed the integration domain from a triangle into a square, see Fig. [2] 
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2.3 A Schwarzschild Black Funnel 



For a Schwarzschild background, consider the following ansatz in d — 5 dimensions: 



ds 2 



xyil + x)' 



x(l-y)(l + xy)Tdt 2 + 



2 x(l + x) 2 Ady 2 



+ 



4y(l-y)(l + xy) 



r 2 B[dx + x(l-x) 2 Fdy} 2 
+ x(l - x) 4 + 



r ' :S -,dn\\, (2.4) 



(l-a;)* 



where T, A, B, F, and S are all functions of x and y. Also, r = 1/2, a choice that we will 
explain shortly. We choose the line element of our reference metric by setting T = A = B = 
S = 1, and F = 0. 

Now let us discuss the boundary conditions. At x — 1, we would like our metric to 
approach that of the planar black hole. This would mean that far from the boundary 
black hole in the field theory, the stress tensor becomes a thermal fluid. Imposing Dirichlet 
boundary conditions, namely T — A — B — S — 1, and F — 0, the line element then reduces 
to: 

ds 2 



y 



dy 2 



+ 



r\dx 2 



■ 4 (1 y 2 )dt 2 + ^ _ ^ , ^_ x y , 4(^^)2 



+ 



(2.5) 



Under the coordinate transformation t = 2r, x = 1 — r /2R and y = z 2 the former line 
element becomes 



-(l-z 4 )rfr 2 + 



dz 2 



+ di? 2 + i? 2 rffi 2 



(2.6) 



which can be recognized as the line element of a five-dimensional planar black hole. 

Our conformal boundary lies at y — 0. If we again choose T — A — B — S — 1, and 
F = 0, the metric becomes 



ds 2 = 



y 



4y x(l + x) 



d^ 2 

2 aS d 



(2.7) 



dsl = -x dt 2 + 



, i 4 + 77-^S^2 , (2-8) 
x(l — x) A (1 — x) 2 

where ds 2 d is the line element of our conformal boundary metric. With the coordinate trans- 
formation x = 1 — r /r, it becomes a more familiar form of the Schwarzschild metric: 



ds 2 9 = 



To 

r 



dr 2 



)dl 2 ■ | "' , ■ rdil . 



(2.9) 
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Of course, the Schwarzschild black hole has the temperature 



(2.10) 



The horizon of our funnel lies at y = 1, where we impose regularity. Note that at y = 1, 
the dt 2 and dy 2 terms of the metric are the same if T = A. This is necessary in order for 
the boundary conditions to be consistent with £ = 0, given our choice of reference metric. 
Expanding the equations of motion about the horizon will give the condition T = A as well 
as other conditions on d y A\ y= i, d y B\ y= i, d y F\ y= i, and d y S\ y= \. 

The temperature associated with the funnel horizon is given by 



T 



1 

2^ 



(2.11) 



We would like the temperature of our black funnel to match the temperature of the boundary 
black hole. In the field theory, this means that the black hole must have the same temperature 
as its surrounding plasma, ensuring that we are constructing the gravitational dual to the 
Hartle-Hawking state. Therefore, (2.10) and (2.11 ) imply that ro = 1/2, justifying the choice 



made below (2.4). 



Finally, there is another boundary condition at x = where we would like our solution 
to approach the zero-energy hyperbolic black hole (2.3). We impose the same Dirichlet 
boundary condition T = A = B = S = 1, and F = 0, and the metric becomes: 



ds 2 



I? 

y 



{l-y)dt 2 



dy 2 



dx 2 1 - 

— ; r H n H dVti 

4y(l - y) Ax 2 4x 



(2.12) 



where we have substituted tq = 1/2. Both Eq. (2.12) and Eq. (2.3) should agree with each 
other close to x — > (rj — > +oo). In order to see this, consider the following change of 
variables y = z 2 and x = exp(— 2rj). Under this coordinate transformation (2.12) reduces to 



ds 2 



^2 



-(l-z 2 )dt 2 + 



dz 2 



l-z< 



+ dr] 2 + 



= 2'/ 



-dm 



(2.13) 



which matches (2.3) as rj — > +oo. 
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2.4 Black Funnels with d = 4 

Now let's consider the d — 4 ansatz: 
L 2 



ds A = 



x 2 y 2 g(x) 2 



Vi \ f/ ^,.2 , x 2 g{x) 2 Ady 2 
x (1 -y)f(x,y)Tdt 



+ 



45(cte + x(l - x 2 ) 2 Fdy) 



;i -z 2 ) 4 



with 



(2.14) 



/(x,y) = 1 + j/ + xV(3-2x) #(x) = 2 + x 2 (3-2a;) £(x) = 1 + A(l - x 2 ) . (2.15) 

As before, T, A, B, F, and 5 are all functions of x and y. The reference metric again has 
T — A — B — S — 1, and F = 0. The constant A controls the size of the boundary black 
hole. The factors of 3 — 2x were added to improve numerics by reducing steep gradients. 

The boundary conditions are essentially the same as the Schwarzschild funnel. At x — 1, 
we require T — A — B — S — 1, and F = and recover the planar black hole: 



ds 2 



1 3 2 dy 2 Adx 2 d<p 2 

--{l-y)dt + r -^ + 9(1 _ x2)4 + 9(1 _ l2)2 



(2.16) 



This can be put in a more familiar form by the transformation t = 3r, x = a/1 — 1/3R, 
which at large R becomes 



ds 2 



-(l-y*)dT 2 + 



dy 2 



+ dR' + R 2 d(f) 2 



l-y3 



(2.17) 



At the conformal boundary y = 0, we again choose T — A — B — S — 1, and F = 0: 



ds 2 



<isg = x 2 <it 2 



1 ds 2 

x z g{x) z 



Adx 2 



l + \(l-x 2 ) 

'i-x 2 y + (i -x 2 ) 2 



d(f) 2 , 



(2.18) 
(2.19) 



With the coordinate transformation x — > y/l — 1/r, our boundary metric becomes 

rfs 2 = - A - ±) rft 2 + ^ + r 2 (l + ^ rf0 2 . (2.20) 
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This is the metric of an asymptotically flat three dimensional black hole with horizon at 
r = 1. The size of the black hole is given by the size of the circle at the horizon: a/1 + A. 
This black hole is not Ricci flat and does not solve the vacuum Einstein equation. Here, 
the boundary black hole is merely a boundary condition and does not need to satisfy any 
equations of motion. In the field theory, this amounts to choosing a fixed background. 

At y = 1, we have a horizon where we require regularity. As before, the dt 2 and dy 2 
terms are the same if T = A. The temperature of this black hole is given by 

T= — , (2.21) 

which is the same temperature as our boundary black hole. 

At x = 0, we again choose the same boundary condition T — A — B — S — 1, and 



F = 0, and the metric becomes a four-dimensional zero energy hyperbolic black hole (2.3) 
after the coordinate redefinition t = 2r, y = z, and x = a/1 + A exp(— 77): 



T 2 

ds 2 = ^ 

z 2 



dz 2 e 2ri 



(2.22) 



2.5 Numerics 

To solve the PDEs numerically, we approximate the system by a set of non-linear algebraic 
equations using pseudospectral collocation on a Chebyshev grid. We then solve the algebraic 
equations with a standard Newton-Raphson relaxation procedure. This method requires a 
seed, which we choose to be the reference metric. 

This method has the advantage that it achieves exponential convergence as the number 
of grid points is increased. We only require a few points (around 40 x 40), and modest 
computational resources. However, a drawback is that the solutions must be analytic. The 
presence of logarithms from a conformal anomaly prevent us from varying the size of the 
black funnel in five dimensions. 

As mentioned earlier in this section, Ricci solitons have been proven not to exist for our 
system. We therefore use the norm of the DeTurck vector \ = to monitor our numerical 
error. In Fig. [3j we plot the maximum value of x as a function of grid points and show the 
expected exponential convergence. 
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Figure 3: The maximum value of \ as a function of the number of grid points N for the 
Schwarszchild funnel (left) and the four-dimensional funnel with A = (right). The 
vertical scale is logarithmic, and the data is well fit by an exponential decay. There is 
some noise in the convergence below 10 -11 from machine precision. 



3 Boundary Stress Tensor 



3.1 Schwarzschild Funnel Stress Tensor 



We now compute the stress tensor for our solutions. We will closely follow de Haro, Skenderis 
and Solodukhin in [12] . We will start by analyzing the five-dimensional funnel. As a first 
step, we need to solve the field equations in a series expansion off the conformal boundary, 
located at y = 0. In order to compute the holographic stress energy tensor in d = 5, we only 



need to solve Eqs. (2.2) to 0(y ). Up to this order, we find the following expansions 



T(x, y) = l + hl- x)x[(l - x)x(3x 2 - 8x + 6) + 1] y 

(1 - x)x[xL T (x) - 4(1 + x)b' 2 (x)} + (6 + Ax - 34x 2 )b 2 (x) 
2(1 -3x)(l + x) 



y 2 + 0(y 3 ), (3.1a) 



A(x, y ) = 1 - 1(1 - x ) 2 xL A (x) y 2 + 0(y 3 ) , 
F(x, y) = (1 - xfx 2 {\ + x)L F (x) y + 0(y 2 ) 



x 



B(x, y) = l- ^(27x 5 - 46x 4 - 23x 3 + 76x 2 - 35x + 1) y + b 2 {x) y 2 + 0{y 3 ) 



(3.1b) 
(3.1c) 
(3.1d) 
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S(x,y) = 1 - -(1 - x)x[(l - x)x(U + 8x- 21x 2 ) + l]y 

4(1 + 2x - 7x 2 )b 2 (x) - (1 - + 4(1 + x)b' 2 (x) - xL s (x)+] 2 3 
+ 4(l-3x)(l + x) V J, ^.lej 

where L^, L^(x), Lp(x) and Lg(x) are polynomials which are given in Appendix [A} This 
expansion is compatible with £ = up to 0(y 3 ). The only free function in this expansion is 
62(2?), which can be extracted by taking two derivatives of B with respect to y. In this way, 
we have no need to resort to a fitting procedure, thus increasing our numerical accuracy in 
determining the holographic stress energy tensor. 

The next step is to change to the familiar Fefferman-Graham coordinates, in which case 
the metric can be recast as 

r2 

ds 2 = — \d~z 2 + dsl + ~z 2 h 2 + z 4 h 4 + 0(z 6 )] , (3.2) 

z 1 

where ds 2 d is defined in (2.9). This coordinate transformation is performed in an expansion 
in 5, where the successive terms are determined by requiring = 0. The coordinate 
transformation has the following schematic form: 

z 2 

+ R 4 (x)2 4 + R 6 (S:)z e + 0(z 8 ) 



x(l + x) 2 

(3.3) 

X = X+ ( 1 -^) 4 ( 1 + 3 ^) g2 + + 0{ ~ 6) 

1 + X 

where the remaining Ri and Vi are given in Appendix |A| The z 2 term in the expansion of y 
fixes the conformal frame we are interested in, namely one where the boundary metric has 
9tt = -x. 

After reading off h^, we can reconstruct the corresponding holographic stress energy 
tensor [12] via 

< T «> - *k ■ (3 - 4 > 

where i and j now run over the coordinates in the boundary metric. Note that because our 
conformal boundary metric is Ricci flat, there is no conformal anomaly and h 2 is identically 
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zero. The holographic stress energy tensor is then given by: 



4nG 5 {Tij) = diag 



x[xf tt (x) - 16(1 - x 2 )b' 2 (x)} + 8 (3 + 2x - 17x 2 ) b 2 (x) 

(3.5) 



8(1 -3x)x(l + x) 5 
x 2 f^(x) - 8b 2 {x) x[xf n (x) - 8(1 - x 2 )b' 2 {x)} + 8(1 + 2x - 7x 2 )b 2 {x 



32£ 3 (l-£ 2 ) 4 ' 32(1 -3x)x 2 (l +x) 3 (l - x 2 ) 2 



This stress energy tensor is manifestly traceless and transverse with respect to (2.9), as 
expected for a conformal field theory on a curved Ricci flat background. are polynomials 
whose expressions are given in Appendix [Aj 

Finally, in Fig. [4] we show how (T t *) varies with r, where r is the radial direction in 



Schwarzschild coordinates defined in Eq. (2.9). The dashed line represents the value of 
(T t *) evaluated in the planar black hole geometry. It is reassuring that the numerical data 
(represented by points) approaches the dashed line for large r. For r/ro > 10, we fit this 
component of the stress energy tensor to a function of the form 

where the last term corresponds to the value of (T t *) for the planar black hole. A x 2 fit yields 
Aq pa —0.007, Bo pa 2.181 and a ~ 0.986, which is consistent with a 1/r tail in the stress 
energy tensor. This should be contrasted with the faster fall off for the same component of 
the stress energy tensor found in [TU]. The slower falloff for the funnels is indicative of the 
stronger coupling between the boundary black hole and its surrounding plasma. Finally, we 
note that (T t *) is smooth across the boundary black hole horizon, as expected. 

3.2 d = 4 Funnel Stress Tensor 

Now we repeat this construction for the four-dimensional funnels. An expansion of the 
equations of motion off the conformal boundary gives 

T(x, y ) = 1 + 1(1- x) 2 K T (x) y 2 + t 3 (x) y 3 + 0(y 4 ) , (3.7a) 

A(x } y) = l + 0(y 4 ), (3.7b) 

F{x, y) = ^(l- x)x 2 (l - x 2 ) 2 g(x)K F (x) y 3 + 0(y 4 ) , (3.7c) 
16 
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Figure 4: Plot of (T t *) as a function of r/ro- The dashed line indicates the value of (T t *) computed 
for the planar black hole, while the red dots represent our numerical data extracted after 



solving (2.2) with appropriate boundary conditions. This particular run was done with 



41 points along each Chebyshev grid. 



B(x, y) = 1 + \{1 - x) 2 K B (x) y 2 + b 3 (x) y 3 + 0(y 4 

o 

1 



S(x,y) 



[l-xyK s (x)y 2 + s 3 (x)y' + 0(y 4 



(3.7d) 



(3.7e) 



where K?, Kp, Kb, and K$ are known functions denned in Appendix A. Furthermore, £3 
and S3 can be expressed in terms of 63 via 



* 3 (x) = -s 3 (x) - 63 (x) , 



S3(X 



63 fx) - 



x(l - x 2 )£(x) [4(1 - x) 3 x(l + 2x) - 18(1 - x)x6 3 (x) + g(x)b' 3 (x)} 
g(x) [x 2 (l - x 2 )X + (1 - 3x 2 )£(x)] 



(3.8a) 
, (3.8b) 



where we have defined £ and g in (2.15). 



Now we go to Feffer man- Graham coordinates: 

L 2 



z 2 



[dz 2 + ds 2 a + z 2 h 2 + z 3 h 3 + 0(z 6 



(3.9) 
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where dsg is given by (2.19). The coordinate transformation that does this is 



xg[x) bx z g{x) q: 



(l-x 2 ) 4 [2 + x 2 (9-8x)] 2 TT . . , , nl fiN 
8xg{x) 



(3.10) 



where Q 3 and C/4 are given in Appendix A. The holographic stress energy tensor is then 
given by [12J 

< T «> " ' < 3 -"> 

Although the four-dimensional funnels are not Ricci flat, there is no conformal anomaly in 
this dimension. Then the boundary stress tensor is 

167rG 4 f 2g(x)-3t 3 (x)-A 4(g(x) + 3b 3 (x) - 2) £ (x)[g(x) + 3s 3 (x) - 2] 

3 { l3 ' iag \ 3xg(x) 3 ' 3x 3 (l-x 4 ) 4 g(x) 3 ' 3x 3 (l - x 2 ) 2 g{xf 

(3.12) 



Combined with (3.8), this stress tensor is manifestly transverse and traceless. 



In Fig. [5j we plot (T t l ) as a function of the radial coordinate r defined in (2.20), for 
various values of A. The behaviour at large r also approaches that of the planar black hole, 
which has the value shown by the dashed line. As before, we can fit this component of the 
stress energy tensor for all points r > 10 to 

We find that for all values of A between and -0.9, A w 0.0021, B w 1.8, and a w 0.99. 
The four- dimensional funnels give the same 1/r falloff as those in five dimensions. This is 
expected since both boundary metrics approach Minkowski space asymptotically as 1/r. 

We also see that (T t l ) is smooth across the boundary black hole horizon. In addition, as 
the size of the black hole is decreased, the value of (T t *) increases at the horizon. Moreover, 
for small black holes, i.e. A ~ —1, (T t *) at the horizon behaves as TZ 3 ^ 2 , where 1Z represents 
some measure of the curvature, such as the square root of the Kretschmann scalar at the 
horizon of the three-dimensional boundary black hole. 
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1 2 3 4 5 



Figure 5: Plot of (T t *) as a function of r for the four- dimensional funnel with A taking values of 0, 
—0.6, —0.8, and —0.9. Lowering A increases (T t *) at the horizon r = 1. The dashed line 
indicates the value of (T f *) computed for the planar black hole, while the coloured dots 
represent our numerical data. We used a 40 x 40 Chebyshev grid. 



4 Embedding into Hyperbolic Space 

In this section, we will visualize the horizon of our funnels by embedding them in hyperbolic 
space. Regardless of the dimension, the induce metric on the spatial horizon of our funnels 
has the same schematic form: 



ds 2 



L 2 



x 



~a (d) 

yxx 



x) dx 1 



X 



(4.1) 



where g x d x(x) and g^'ix) are functions that depend only on x and can be read from the 



pullback of the line elements (2.4) and (2.14) to the y = 1 hyperslice. 



On the other hand, d — 1 Euclidean hyperbolic space has the following form 



ds 



-^{dy 1 + dR 2 + R'dnj„ 3 ) . 



(4.2) 



An isometric embedding can be easily achieved if we consider a curve of the form s(x) 



(y(x), R(x)) in hyperbolic space. The pullback of the line element (4.2) to s gives the 
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following d — 2 induced metric: 



Demanding that (4.3) agrees with (4.1) gives a systems of two equations in R(x) and 
y(x). One of the equations is a simple algebraic equation that determines R(x) as a function 
of y(x), whereas the remaining equation is a first order non-linear ODE in y{x) which we 
solve numerically on a Chebyshev grid. As a boundary condition, we fix -R(O) to be the size 
of the boundary horizon in the black hole frame. Note that hyperbolic space is appropriate 
for the visualization of our funnels since a Schwarzschild string in AdS would be represented 
by a line of constant R, and the planar black hole would be given by a line of constant y. 

The embeddings are displayed in Fig. [6) The left panel shows the Schwarzschild funnel 



while the right panel shows the four-dimensional ones. The vertical dashed line in Fig. 6(a) 
represents an embedding digram for a Schwarzschild string with r = 1/2. One can attempt 
to give an heuristic argument based on these embeddings in favor of the stability of the 
Schwarzschild funnel. The stability of the Schwarzschild string in AdS was first studied in 
[T5] where it was found that the string is Gregory-Laflamme unstable only when it exceeds 
a length of about Ay rj 2.94. One can ask when our geometry significantly deviates from 
the AdS-Schwarzschild string} In order to quantify this, we define AR = |1 — R(x)/tq\. 
We do not consider the black funnel to be string-like if we find deviations above the five 
percent level, or equivalently when AR < 0.05. Given that this occurs at about y ~ 0.53, 
we find it unlikely that this geometry will be Gregory-Laflamme unstable, since the unstable 
mode does not seem to fit the string-like region of the Schwarzschild funneQ Indeed, we 
have attempted to directly study the stability of these objects to axisymmetric perturbations 
and our preliminary results indicate stability in this sector of the perturbations. We have 



also searched for negative modes and found none. In Fig. 6(b) one sees that the four- 
dimensional case is even less likely to be unstable, given that the four-dimensional funnels 
appear to widen more quickly than the five-dimensional case. Even if the size is decreased, 
the four-dimensional funnels do not become more narrow. Furthermore, there is no known 
Gregory-Faflamme instability in four dimensions. 



3 Note that a string is unstable if the mode that sits at the onset of the Gregory-Laflamme instability fits 
inside the string. 

We thank Donald Marolf for this heuristic argument. 
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R(x) 

(a) d = 5 



R(x) 

(b) d = 4 



Figure 6: Embedding diagrams of the spatial horizon induced metric for d = 4 (left) and d = 5 
(right). The dashed line represents the embedding of a Schwarzschild black string horizon. 
The values for A in the d = 4 plot are 0, -0.6, -0.8, -0.9, and -0.99. 

5 Discussion 

We have numerically constructed black funnel phases dual to strongly coupled Hartle- 
Hawking states on a Schwarzschild background and on a class of three dimensional black hole 
backgrounds. We computed their holographic stress energy tensors and also varied the size 
of the three dimensional black hole backgrounds. We also embedded the horizon geometry of 
these funnels in hyperbolic space, and argued that no Gregory-Laflamme instability occurs 
for these solutions. 

There are a number of future directions. First, it would be interesting to complete the 
stability analysis of the black funnels, i. e. consider non-axisymmetric perturbations. Even 
though we don't envisage any other type of instability, a rigorous analysis would be most 
welcome. 

Here, we have only varied the size of the black hole in four dimensions where there is 
no known Gregory-Laflamme instability. It would be interesting to see if such an instability 
occurs in higher dimensions. One could attempt to find six-dimensional funnels where there 
is no conformal anomaly, or to vary the size of five-dimensional funnels using some other 
numerical method like finite differences. 

In order to fully test the phase transition conjectured in [5], the black droplet phases 
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must also be constructed. These solutions naturally lie in a domain of integration with five 
boundaries: two horizons, an axis, the planar black hole metric, and the boundary metric. 
This type of problem likely requires patching of separate integration domains, or some means 
of collapsing one of these boundaries to a point. Three of these boundaries are also 'fictitious' 
and require Neumann data rather than Dirichlet, making it more difficult to find a suitable 
seed. Perhaps one can attempt to find solutions sufficiently close to the droplet in [10], or 
obtain a suitable seed from a matched asymptotic expansion. 

Our solutions have bulk horizons that are at the same temperature as the boundary 
horizons. One can imagine the field theory scenario where the boundary black hole is kept 
at a different temperature than its surrounding plasma. In order to construct the bulk dual 



to this situation, one can use a different hyperobolic black hole than (2.3). 

One could attempt to find equivalent solutions in global AdS, as opposed to our funnels 
which can be thought of as lying in the Poincare patch. If the boundary field theory lies on 
a sphere with two black holes at the poles, there are now potentially three competing phases 
for the gravitational dual. There is a funnel phase where the two boundary horizons are 
connected in the bulk. But there are also two droplet phases where the boundary horizons 
are disconnected: one with a third spheroidal horizon in the bulk, and one without. 

Aside from Hawking radiation, black funnels have also been used to study heat transport 
[!%] . A black hole is viewed as a heat source in the dual field theory and the heat transport 
properties are mapped to a stationary flow in the black funnel. One might be able to 
construct stationary black funnel solutions in order to study non-equilibrium heat flow. 
This can also be attempted in global AdS. 

One can also imagine having a boundary black hole that rotates such as the Kerr metric. 
Then the dual black funnels will also be rotating. In order to prevent superluminal motion 
far from the axis of rotation, the horizons must 'twist'. These solutions will therefore have 
an event horizon that is not also a Killing horizon^} 
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A Asymptotic expansion 

In this appendix we present the functions used in section 3. The functions for the Schwarzschild 



funnel appearing in (3.1), (3.3) and (3.5) are: 



L T (x) = 4185a; 11 - 8409a; 10 - 13770a; 9 + 44406a; 8 - 16301a; 7 - 40548a; 6 + 37538a; 5 

+ 2401a; 4 - 15024a; 3 + 6411a; 2 - 913a; + 26 , (A.la) 

L A (x) = -81a; 9 +418a; 8 -1037x 7 + 1336x 6 -466a; 5 -840a; 4 + 1002x 3 -356x 2 + 23x + 2, (A.lb) 

L F (x) = -27a; 4 + 31a; 3 + 22a; 2 - 30a; + 5 , (A.lc) 

L s (x) = 5400a; 11 - 11460a; 10 - 14496a; 9 + 51308a; 8 - 18576a; 7 - 46818a; 6 + 41490a; 5 

+ 4333a; 4 - 17163a; 3 + 6892a; 2 - 939a; + 32 , (A. Id) 

„ , , 95 6 - 305 5 + 315 4 - Ax 3 - 95 2 + x + 2 . , 

RAx) = j , (A.le) 

V ; 25 2 (5 + l) 4 V 1 

Rq{x) = [8b 2 (x) - 5 2 (23495 10 - 106445 9 + 151505 s + 644a? 7 - 218175 6 + 181465 5 

+ 13425 4 - 91765 3 + 50235 2 - 11145 + 95)]/[325 3 (l - 5 2 ) 4 ] , (A.lf) 

- . (1 - 5) 4 (725 6 - 645 5 - 1765 4 + 1925 3 + 405 2 -635-3) , . , 

V 'V = " ~ 4(1+1? 1 ' (A ' lg) 

f tt = -164435 12 + 485225 11 + 251675 10 - 2331805 9 + 2344455 s + 1077325 7 

- 3132035 6 + 1326905 5 + 756375 4 - 865325 s + 285785 2 - 35325 + 95 , (A. lh) 



fss = 23495 10 - 106445 9 + 151505 8 + 6445 7 - 218175 6 + 181465 5 

+ 13425 4 - 91765 3 + 50235 2 - 11145 + 95 , (A.li) 
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f n = -11745x 12 + 378785 11 + 1677x 10 - 138028x 9 + 1568795 s + 48786x 7 



187669x e + 87840x 5 + 40131x 4 - 51206x 3 + 177725 2 - 2418x + 95 . (A.lj) 



^3 



Now for the four-dimensional funnel, we give the expressions for the functions appearing 



in (3.7), (3.10). They are given by 



K T {x) = 4[l + x(2 + 9x 3 (l 



-2\4 



X 



+ {l-x 2 ){l + x) 2 H^{x) 



(A.2a) 



K F (x) = 24(1 - xf{\ + 2x) + 3(1 - x 2 fH + {x) - 6(1 - x)x\l + x) 2 g(x^ 



(1 - x)x 2 (l + x) 2 (l - x 2 )X 2 J(x)g(x) 
J(x) 2 



+ 



3x 2 (l — x) 2 pi(x) 
I(x) 



+ 3p 2 {x) 



g(x) (A.2b) 



Q 3 



U A {x) 
where 



K B {x) = Ap 3 {x) - (1 - x 2 ){l + x) 2 H_(x) (A.2c) 

K s {x) = 12(1 - x 2 ) V [6(1 - x 2 ) 2 x 2 + (1 + a;) J(a;)] - K T {x) (A.2d) 

4(1 - x) 2 [l + x(2 + 9x 3 (l - x 2 ) A )] + (1 - x 2 ) 4 c/(x) [12(1 - x)x 2 + #(»] 2 ^ 

16a; 3 5'(x) 3 

4p 4 (x) + (l-x)[6(l-x)x 2 + ^(x)]ir i j(x) + 2x 2 J K" F (x) [ , (A.2f) 



(1 -x)(l -x 2 ) 
128x 3 c/(x) ; 



i? ± (x) = 
J(x) 



^(x) 2 
£(x) 2 

#(x) 



x 4 (l-x 2 )A 2 ±(l + 3x 2 )£(x) 2 



(A.3a) 



(x) L 



P2(^) 



x 2 (l-x 2 )A + (l-3x 2 )£(x) (A.3b) 

pi(x) = 32x 5 - 43x 4 - 17x 3 - 8x 2 - 2x + 8 (A.3c) 
18x n -57x 10 +x 9 -38x 8 -112x 7 +327x 6 + 157x 5 -316x 4 -72x 3 +88x 2 +8x-2 (A.3d) 



p 3 (x) = -33x 12 + 15a; 11 + 126x 1 + 15x 9 - 114x 8 - 81x 7 - 30x 6 + 57x 5 + 63x 4 

- 6x 3 - 12x 2 - 2a; - 1 (A.3e) 

p 4 (x) = 6(l-x)V(l+2a;) + {l+x(l+x) l+(l-x) 2 x(l+x) 2 (2+x 2 (9-8x)) 2 }#(x) , (A.3f) 



and (7 and £ were already defined in (2.15) 
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